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We investigate pairing symmetry in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) vortex state 
and vortex lattice, and explain the electronic structure in these states in terms of pairing symmetry. 
We show analytically that at the intersection point of FFLO nodal plane and vortex line, only even 
frequency pairing is present if the Zeeman splitting is negligibly small. With increasing Zeeman 
splitting, odd frequency pairing emerges there. This makes it possible to interpret the gap structure 
of the density of states at the intersection point as a manifestation of the even frequency pairing. 
In the vortex lattice, we find that only odd frequency pairing is present at the core centers, while 
at the midpoint of the vortex lines, only even frequency pairing exists. Thus, the odd and even 
frequency pairings also form the lattice in the vortex lattice state. 



I. INTRODUCTION 

The study of the mixed state in type-II supercon- 
ductors has a long history and revealed a wide variety 
of physical phenomena.^ In the clean limit, low-energy 
bound states, dubbed Andreev bound states, are formed 
in the vortex core due to the spatial profile of the su- 
perconducting pair potential. The bound states mani- 
fest themselves as an enhancement of zero-energy quasi- 
particle density of states (DOS) in the core, observable 
as a zero-bias conductance peak by scanning tunneling 
microscope. Recently, the Andreev bound states have 
been revisited from the viewpoint of the odd frequency 
pairing.^ 

Generally, superconducting pairing is classified into 
even-frequency or odd-frequency state according to a 
symmetry with respect to time. Due to the Fermi statis- 
tics, even-frequency superconducntors belong to the sym- 
metry class of spin-singlet even-parity or spin-triplet odd- 
parity pairing state, while odd- frequency superconduc- 
tors belong to the spin-singlet odd-parity or spin-triplet 
even-parity pairing state. 

Although the possibility of the odd-frequency pair- 
ing state in various kinds of uniform systems was dis- 
cussed in the literature its realization in bulk materi- 
als is still controversial. On the other hand, the realiza- 
tion of the odd-frequency pairing state has recently been 
pointed out in inhomogeneous even-frequency super- 
conducting systems such as ferromagnet /superconductor 
junctions,^ normal metal/superconductor interface^ii^ 
Abrikosov vortex^ or a diffusive normal metal attached 
to a spin-triplet superconductor.^^ 

In Ref. 5, pairing symmetry in an isolated vortex is elu- 
cidated. It is found that at the center of the core, purely 
odd frequency pairing state exists. Since single vortex is 
considered in Ref. [Bl, this theory is applicable to the low 
field regime. On the other hand, at high magnetic fields, 
the distance between vortices becomes short, and hence 
the overlap effect of the vortex cores cannot be neglected. 
In this situation, the local DOS around a vortex core is 



expected to break cylindrical symmetry and show six- 
fold symmetric structure when the vortex lattice forms a 
triangular lattice. ^^^^^ Also, at high field regime, Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) vortex state may be 
realized under some conditions. These states are also 
inhomogeneous systems, which may provide rich struc- 
ture of superconducting correlation, in particular odd 
frequency pairing. 

The FFLO state^i^ is realized by Cooper pairs with 
non-zero center of mass momentum when the Fermi sur- 
faces for up-spin and down-spin electron bands are largely 
split by the Zeeman effect, thus inducing the spatial mod- 
ulation of the pair ^oieniml}^ '^^'^•^^•^^•^ Several 
experiments support the realization of the FFLO state 
in a high field phase of a quasi-two dimensional (Q2D) 
heavy Fermion superconductor CeCoIns .^^ ^^^^^^^^^^^^i^^i^-^ 
There, it is supposed that the FFLO nodal planes of the 
pair potential run perpendicular to the vortex lines. 

In general, we have to include the vortex structure in 
addition to the FFLO modulation, because the FFLO 
state appears at high fields in the mixed states i^-'^i^^i^^ 
The vortex and FFLO nodal plane structures in the 
FFLO state were studied in Refs. [32II33I . It is found that 
the topological structure of the pair potential plays im- 
portant roles to determine the electronic structures in 
the FFLO vortex state. The pair potential has 27r-phase 
winding around the vortex line, and 7r-phase shift at the 
nodal plane of the FFLO modulation. These topologies 
of the pair potential structure affect the distribution of 
paramagnetic moment and low energy electronic states 
inside the superconducting gap. 

Another aspect of the FFLO state is the parity mixing. 
In the FFLO state, in addition to the breakdown of the 
translational symmetry, Zeeman field breaks the SU(2) 
symmetry in spin space. As a consequence, singlet-triplet 
mixing state emerges as a stable phase in the FFLO 
stoXeMi^^^ Therefore, quite gerenally, in the FFLO 
state mixture of even-odd frequency and singlet-triplet 
pairings is expected to emerge. 

In this paper, based on the quasiclassical theory of su- 
perconductivity, we investigate pairing symmetry in the 
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FFLO vortex and vortex lattice, and explain the elec- 
tronic structure in these states in terms of pairing sym- 
metry. We show analytically that at the intersection 
point of FFLO nodal plane and vortex line, only even 
frequency pairing is present if the Zeeman splitting is 
negligibly small. With increasing Zeeman splitting, odd 
frequency pairing emerges there. This makes it possi- 
ble to interpret the gap structure of the density of states 
at the intersection point as a manifestation of the even 
frequency pairing. In the vortex lattice, at the core cen- 
ters, only odd frequency pairing is present while at the 
midpoint of the vortex lines, only even frequency pairing 
appears irrespective of the energy. Therefore, the odd 
and even frequency pairings also form the lattice in the 
vortex lattice state. 

The organization of the paper is as follows. In Sec. II, 
we study the FFLO vortex state, and show the electric 
structure and superconducting correlation of the system. 
Sec. Ill is devoted to the study of the vortex lattice. The 
summary is given in Sec. IV. 



II. FFLO VORTEX 



states in this paper, since some experiment al^i^i^ and 
theoretical^^'^^ works support the LO state (at least in 
low temperature region) for the FFLO states in CeCoIns. 

When we consider vortex structure in the LO state, 
there are two possible choices of the configuration for 
the vortex lines and the FFLO modulation: the modu- 
lation vector of the FFLO state is parallel^^ or perpen- 
dicular^^'^^ to the applied magnetic field. In this paper, 
we study the former case by the quasiclassical theory as 
shown in Fig. ^d^^^^dl 
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A. Formulation 

The electronic structure of the vortex core in an in- 
homogeneous clean superconductor can be described by 
the quasiclassical Eilenberger equations^i^ based on the 
Riccati parametrization^^ . Along a trajectory r{x') = 
ro-\-x' vf with unit vector vp parallel to Fermi velocity 
vf^ the Eilenberger equations are generally represented 
in 4x4 matrix form^^, which, in terms of the Ricatti pa- 
rameters a and 6, reduces to: — 

vpdx'CL + 2ujnCL + ttA'^'a — A + iSa — = (1) 
vpO^'h - 2ujnb - bAb + A"^ - ibJ: ^itb = (2) 

with the self energy of the Zeeman splitting S = S = 
IIbHcJz and Pauli matrix (Ji{i = x^y^z). Here, we have 
defined 2x2 matrix d and b in spin space via the 4x4 
Green's functions g: 
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There are two possible spatial modulation of the pair 
potential A in the FFLO states. One is the Fulde-Ferrell 
(FF) state^^ with phase modulation such as A ex e^^^, 
where Q is the modulation vector of the FFLO states. 
The other is the Larkin-Ovchinnikov (LO) state^^ with 
the amplitude modulation such as A ex sin Qz, where the 
the pair potential shows the periodic sign change, and 
A = at nodal planes. We discuss the case of the LO 



FIG. 1: (color online) Configurations of the vortex line and 
the FFLO nodal planes in the x — z plane. 

As a model of Fermi surface in CeCoIns, we use 
a Q2D Fermi surface with rippled cylinder-shape, and 
the Fermi velocity is given by vp = {va^Vb^Vc) cx 
(cos^, sin^, 'D;^ sin /cc) at the Fermi surface kp = 
{ka^ki)^kc) cx (/cpo cos6>, /cpo sin^, /cc) with < 6> < 27r 
and — TT < < TT.^^ In our calculation we set i}^ = 0.5, so 
that the anisotropy ratio 7 = ^c/^ab ^ (^?)k^^/(^a)k^^ ^ 
0.5. Here, vf = (v^)^^'^ is an averaged Fermi velocity on 
the Fermi surface. (•••)!< indicates the Fermi surface av- 
erage. A magnetic field is applied along the a axis direc- 
tion in our calculation. Thus, the coordinate (x, z) for 
the vortex structure corresponds to (6, c, a) of the crystal 
coordinate. 

We consider (i2;2_^2-wave FFLO vortex state, and 
choose the following form of the pair potential: 



A(r,6>) = Aocos26>tanh 



exp(z(/:)) sin QzcFyiA) 



with exp(z(/p) = (x + iy)/ ^x'^ + y'^. Here, we introduce 
the coherence length ^ = f^'Up/Ao, the vortex line is situ- 
ated at X = 7/ = 0, and exp(z(/?) is the phase factor which 
originates from the vortex. 

For the calculation of the local DOS normalized by its 
value in the normal state, the quasiclassical propagator 
has to be integrated over and kz which define the di- 
rection of the Fermi velocity. The normalized local DOS 
in terms of functions d and b is given by 
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iV(ro,^) = -Tr 



— / —Re 



(l + a6)-^(l-a6) 



(5) 



r 



where £^ denotes the quasiparticle energy with respect to 
the Fermi level and S is an effective scattering parameter. 
In numerical calculations throughout this paper, we will 
fix this value as (5 = O.OSAq. 

We obtain singlet (triplet) pair amplitude /s(^) as 



fsm = {f''-{+)f'''}/2. 

Its even (odd) frequency component, 
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is obtained as 



even{odd) 



(6) 
and 



fLenioM) = {fs{0, k.) + HfsiO + 7T, -k,)} /2, (7) 
feveniodd) = iMO, " i + )ft{0 + TT, -k,)} /2. (8) 

Note that due to the Fermi statistics, for singlet pairing, 
even (odd) parity state should be even (odd) frequency 
pairing, while for triplet pairing, even (odd) parity state 
should be odd (even) frequency pairing. 

Their average in momentum space is defined as 
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even{odd) 



1 p'^TT 



'A 

even{odd) 



dOdk^ (9) 



with A = s^t. The magnitude of the even (odd) fre- 
quency pairing can be repesented as 

(feveniodd)) ^ (^flven{odd)) ^ (^flven{odd)) ' i^^) 



B. Results 

In the FFLO vortex state, along the trajectory through 
the intersection point of a vortex and a nodal plane, the 
gap function does not change sign, because the phase 
shift is 27r by summing tt due to the vortex and tt due to 
the nodal plane. ^^^^^ Consequently, the DOS at the in- 
tersection point does not have the Andreev bound state 
which is seen in the conventional vortex state or FFLO 
nodal plane. Our analysis presented below makes it pos- 
sible to understand these features of DOS from the view- 
point of the symmetry of superconducting correlations. 

First, we discuss the general property of the symmetry 
in the FFLO vortex state. Consider a trajectory passing 
through the intersection point of vortex line and FFLO 
nodal plane. By setting = at the intersection point, 
we obtain h{x' ^uOn) = (Jyd{—x'^ —U0n)~^(yy from the Eilen- 
berger equations since the gap function does not change 
sign at = and K{x') = A{—x'). 

Thus, at the intersection point x' = 0, we have 
/(0,cJn) = cryf{0,-ujn)cry, and hence fs{0,ujn) = 
fsiO.-LJn) and /t(0,cjn) = -ft{0,-iJn)' This means 



that at the intersection point, only even- frequency spin- 
singlet even-parity and odd-frequency spin-triplet even- 
parity parings are allowed to exist [Note that even (odd) 
frequency singlet (triplet) pairing should have even par- 
ity in accordance with the Fermi statistics]. In particular, 
at jiisH = since spin is conserved, at the intersection 
point only even- frequency spin-singlet even-parity pair- 
ing exists. With increasing fisH, odd- frequency spin- 
triplet even-parity pairing emerges there. In a similar 
way, we can show that at the center of the FFLO nodal 
plane (vortex core) without vortex (the FFLO modu- 
lation), /{O^cOn) = —(Jyf{0^—u;n)o'y is Satisfied since 
the gap function changes sign at = and A{x') = 
—A{—x^). Hence, only odd-frequency spin-singlet odd- 
parity and even-frequency spin-triplet odd-parity pair- 
ings are allowed there. In particular, at /JjbH=0, only 
odd-frequency spin-singlet odd-parity pairing is present. 
The suppression (enhancement) of the local DOS at zero 
energy is related to the presence of even (odd)-frequency 
pairing. ^^^^^^ It has been shown that the emergence of the 
odd-frequency pairing is a physical reason of zero energy 
peak of the local DOS inside the vortex core.^ 

Next, we show the results of DOS and pair amplitudes 
by numerically solving the Eilenberger equation. In the 
following we set Q = 27t/ (100^) and y = 0. 

Figure [2] shows (a) the local DOS, (b) magnitude of 
even frequency pairing and (c) magnitude of odd fre- 
quency pairing as a function of the energy at the intersec- 
tion point (x, z) = (0, 0). At jiisH = 0, local DOS shows 
a gap structure, which reflects that at the intersection 
point, only even frequency pairing exists at /i^^ = as 
seen in Figs. [2](b) and (c). With increasing fisH, these 
structures are split and odd frequency pairing emerges 
there. As shown in the above discussion, at the intersec- 
tion point, only even parity pairing is allowed to exist. 
Thus, even (odd) frequency pairing should be spin singlet 
(triplet). 

Figure [3] displays (a) the local DOS, (b) magnitude 
of even frequency pairing and (c) magnitude of odd fre- 
quency pairing as a function of the energy at the center 
of the core, {x,z) = (0,250- At /j^bH = 0, local DOS 
shows a peak structure, reflecting the presence of the odd 
frequency pairing for /j^bH = as seen in Figs. [3](b) and 
(c). As fiBH increases, these structures are shifted and 
even frequency pairing emerges there. Note that here 
spin singlet (triplet) state dominates spin triplet (sin- 
glet) state in even (odd) frequency pairing. This can be 
understood by considering the vortex without the FFLO 
nodal plane. In this case, only odd parity pairing is al- 
lowed at the core center. Therefore, at the core center 
far away from the FFLO nodal plane, we can also expect 
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FIG. 2: (color online) Normalized local DOS at the intersec- 
tion point (x^z) — (0,0). 



FIG. 3: (color online) Normalized local DOS at the the center 
of the core (x, z) — (0, 25^). 



that odd parity pairing dominates. 

Figure m depicts the results at the FFLO nodal plane, 
(x,z) = (25^,0). The results are qualitatively similar to 
those in Fig. [3l Quantitative difference may originate 
from the dimensionality of the Andreev bound states: 
ID in vortex line and 2D in the FFLO nodal plane. The 
local DOS at iibH = shows a zero energy peak by 
the presence of the odd frequency pairing for /j^bH = 
as seen in Figs. 3] (b) and (c). Finite (but small) 
magnitude of the even frequency pairing at jiisH = 
in Fig. m (b) is physically due to the penetration of the 
even frequency pairing at the intersection point into the 
FFLO nodal plane. This results from the fact that the 
relation A{x') = —A{—x') does not hold any more at the 
the FFLO nodal plane due to the presence of the vortex. 
With increasing jubH^ these structures are split and even 
frequency pairing is enhanced. 

Figure [5] shows the results far away from the nodal 
plane and the core (x, z) = (25^,25^). The results are 
qualitatively similar to those in Fig. [2l At jj^bH = 0, 
local DOS shows a gap structure with dominant even 
frequency component for /i^i^ = as seen in Figs. [5](b) 
and (c). With the increase of fiBH, these structures are 
split and odd frequency pairing is enhanced. 

In Fig. [6l we show (a, d, g) normalized local DOS, (b, 
e, h) magnitude of even frequency pairing, and (c, f, i 



) magnitude of odd frequency pairing at /j^bH = with 
(a-c) E = 0, (d-f) E = 0.5 Ao and (g-i) ^ = Aq. At 
E = 0, the DOS shows a strong peak at the FFLO nodal 
plane and vortex line but at the intersection point, it 
is suppressed. Correspondingly, large magnitude of the 
odd frequency (spin singlet odd parity) pairing is seen 
at the FFLO nodal plane and vortex line, while at the 
intersection point or far away from both the FFLO nodal 
plane and vortex line, only even frequency (spin singlet 
even parity) pairing exists. With increasing E^ the DOS 
is reduced at the FFLO nodal plane and vortex line while 
it increases away from the FFLO nodal plane and vortex 
line. The structures of even and odd frequency pairings 
are smeared with the increase of E. 

The corresponding results at /j^bH = O.SAq is shown in 
Fig. [71 It is seen that even and odd frequency pairings 
are spatially distributed in a similar way. At E = 0, 
the DOS has a similar structure to that in Fig. [6] (d) 
due to the Zeeman spliting fiBH = O.SAq. As shown 
in Fig. [71 (b), the structure of even frequency pairing is 
also simlar to that in Fig. [6l (e). On the other hand, 
the behavior of odd frequency pairing in Fig. [71 (c) is 
rather different from that in Fig. [6l (f). This is because 
by introducing the Zeeman spliting, spin triplet pairing, 
especially odd frequency spin triplet even parity pairing, 
emerges, which will be discussed together with Fig. [51 
below. At £^ = 0.5Ao, the DOS has a similar structure 
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large magnitude of even frequency pairing at vortex line 
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FIG. 4: (color online) Normalized local DOS at the nodal 
plane [x,z) = (25^,0). 

to that in Fig. [6] (a). As shown in Fig. [3(f), the structure 
of odd frequency pairing is also simlar to that in Fig. [6] 
(c), while the even frequency pairing in Fig. [71 (e) shows 
a different feature from that in Fig. [6] (b). This also 
results from the generation of the even frequency spin 
triplet odd parity pairing. At ^ = Aq, the results are 
similar to those in Figs. [71 (a-c). 

Now, let us investigate spin and parity structure of 
superconducting correlation. We show (a) magnitude of 
even frequency spin singlet even parity pairing, (b) mag- 
nitude of even frequency spin triplet odd parity pairing, 
(c) magnitude of odd frequency spin singlet odd parity 
pairing, and (d) magnitude of odd frequency spin triplet 
even parity pairing, at £^ = and iibH = 0.5Ao in Fig. 
[HI We see that even parity pairing, {feven) if odd) ^ 
is dominant at the intersection point or far away from 
both vortex line and FFLO nodal plane, while odd par- 
ity pairing, (/g^en) if odds) ^ dominates at vortex line 
and FFLO nodal plane. Therefore, we find that large 
magnitude of odd frequency pairing away from vortex 
line and FFLO nodal plane in Fig. [7] (c) is due to the 
spin triplet pairing induced by the Zeeman term. The 



FIG. 5: (color online) Normalized local DOS away from the 
nodal plane and the core {x,z) — (25^,25^). 



III. VORTEX LATTICE 
A. Formulation 

Here, we will consider vortex lattice in 5- wave super- 
conductor. Following Ref. il2j, let us explain the Eilen- 
berger theory of vortex lattice (see also Refs. I45ll46ll47l ). 
First, we obtain the pair potential self-consistently by 
solving the Eilenberger equation in the Matsubara fre- 
quency. Next, using them, we calculate the local DOS 
by solving the Eilenberger equation in the real energy 
instead of the Matsubara frequency. 

In our calculation, the unit vectors of the vortex lattice 
are given by ri = (aa;,0), Y2 = {(ax, ay). As we consider 
a triangular lattice, we set ay /ax = a/3/2 and ( = 1/2. 
The microscopic magnetic field H = (0, 0, ^) is given by 
H(r) = V X A(r) = H where the vector potential A(r) 
is A(r) = |H X r in the symmetric gauge. 

For the quasi-classical Green functions, the Eilenberger 
equation is given as 
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FIG. 6: (color online) (a, d, g) normalized local DOS. (b, e, h) magnitude of even frequency pairing, (c, f, i) magnitude of odd 
frequency pairing, (a-c) E = 0. (d-f) E = 0.5Ao. (g-i) ^ = Aq. Here, we set /j.bH = 0. 



(^11 + z27rA|| /(l)o)a{ujn,0, r) + {2uJn + A*a(cc;n, 6>, r)) a{ujn, 6>, r) - A = 0, (11) 

(9,1 - i27rA|| /0o)K^n, 0, r) - {2iJn + A^^n, 0, r)) 6(cjn, 6>, r) + A* = 0, (12) 

I 

with the Matsubara frequency cOn = (2n + l)7rT. Here, r is the center of mass coordinate of a Cooper pair. The 
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FIG. 7: (color online) (a, d, g) normalized local DOS. (b, e, h) magnitude of even frequency pairing, (c, f, i) magnitude of odd 
frequency pairing, (a-c) E = 0. (d-f) E = 0.5Ao. (g-i) ^ = Aq. Here, we set /j.bH = O.SAq. 



direction of the relative momentum of the Cooper pair, 
k = k/|k|, is denoted by an angle measured from the x 
axis in the hexagonal plane. Here, we define d\\ = d/dr\\ 

and A\\ = k • A = —^Hr_\_. Also, we have taken the 
coordinate system: u = cos^x + sin^y, v = — sin^x + 
cos , thus a point r = xx + ?/y is denoted as r = ry u + 
r±v. The first-order differential equations (pTj) and (p!2|) 
are solved along the trajectory where r± is held constant. 
Notice that we here focus on superconductor in the type 



H limit, and a weak magnetic field regime so that the 
Zeeman term can be neglected. 

The self-consistent equation for the pair potential A(r) 
reads 

n27r in 

A(r) = yiVo27rT V / — /K,0,r), (13) 

with / = — 2m/(l + a6), the density of states at the 
Fermi surface Aq, and the pairing interaction V. In our 



calculation, we use the relation 
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FIG. 8: (color online) (a) magnitude of even frequency spin 
singlet even parity pairing, (b) magnitude of even frequency 
spin triplet odd parity pairing, (c) magnitude of odd fre- 
quency spin singlet odd parity pairing, (d) magnitude of odd 
frequency spin triplet even parity pairing. Here, we set ^ = 
and /2bH = 0.5Ao. 



We calculate the r.h.s. of Eq. (p!3|) using the quasi- 
classical Green functions obtained by Eqs. (pTj) and (p!2]) , 
and obtain the new value for A(r). Using the renewed 
pair potential, we solve the Eilenberger equation (pTj) and 
(p!2|) again. Using the following gap function as an initial 
value,^- 



A(r) 



1/4 oo 

p= — oo 



ay fy^yo ^ > 

-TT— hp 

ax \ ay ) 



■2m 



'^0 , C 



yo_ 

a,. 



exp I m- 



xy 

axay ^ 



(15) 



we repeat this simple iteration procedure more than 20 
times, and obtain a sufficiently self-consistent solution 
for A(r). In Eq. ([15]), the r.h.s. is the Abrikosov so- 
lution of the vortex lattice, where we use the relation 
Haxay/(j)o = 1. The factor ex.-p{i7rxy /axay) is due to 
the gauge transformation from the Landau gauge to the 
symmetric gauge. We set tq = (xo,^o) = +1*2) 
so that one of the vortex centers locates at the origin of 
the coordinate. Note that since we determine the gap 
function self-consistently, our theory is applicable to any 



magnetic field below the upper critical field, although we 
adopt the Abrikosov solution as an initial value. 

The physical quantities are calculated in a similar way 
to the previous section. The DOS is calculated as 

N{v,E)= / Re[(l + a5)-i(l-a5)]^^^^^^,.(16) 

J 

Even (odd) frequency component feven{odd) is given by 

f even{odd) 

= {}{e) + {-)}{e + n)}/2. (17) 
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Its average is defined as 

2 1 2 

(feveniodd)) ^ 2^ J \feveniodd)\ dO . (18) 

Note that since spin is conserved, all the above pairings 
are spin singlet. 

In the following, we introduce i?o, the transition tem- 
perature Tc, and Hq as units of length, temperature, and 
magnetic field, respectively, where Rq = TiVY/^iik^Tc and 
Hq = hc/2\e\R^. Also, we fix (5 as (5 = O.OSAq where Aq 
is the bulk value of the gap function at T = 0. 

B. Results 

Here, we consider two cases: low field case where the 
distance between vortices is large, and the overlap effect 
of the vortex cores is weak, and high field case as an oppo- 
site situation. In the former case, we choose H/Hq = 0.05 
and in the latter we set H/Hq = 0.5. The temperature 
of the system is fixed as T/Tc = 0.5. 

Figure [9] shows the results at H/Hq = 0.05. At £: = 0, 
the Andreev bound states are seen. The odd frequency 
pairing shows a similar structure to the DOS, while even 
frequency pairing is absent at the core center although it 
has a large value near the core center. As E increases, 
the DOS and the odd frequency component at the core 
center are reduced but they increase away from the core 
center. We find that at the core center, only odd fre- 
quency pairing is present while at the midpoint of the 
vortex lines, only even frequency pairing exists irrespec- 
tive of the energy. Thus, the odd and even frequency 
pairings also form the lattice in the vortex lattice state. 
In Ref. I5I, it is clarified that only odd frequency pair- 
ing is present at the core center of an isolated vortex. 
We see that this is also the case for the vortex lattice. 
We have also found that the orbital pairing symmetry 
of superconducting correlation contains not only p-wave 
component but also higher harmonics, such as d-wave or 
/-wave components, in the vortex lattice since the ro- 
tational symmetry is broken, in contrast to the single 
vortex case.^ 

We show the results at H/Hq = 0.5 in Fig. [TOl A 
qualitatively similar tendency to Fig. [9] is seen. However, 
compared to Fig. [9l the vortex spacing is reduced and 
the overlap effect of the vortex cores reduces the Andreev 



bound states. Also, the magnitudes of the even and odd 
frequency components are suppressed. Meanwhile, DOS 
away from the core increases, and the mangitude of the 
odd frequency pairing becomes comparable to that of the 
even frequency pairing. 

To see how even and odd frequency pairings are influ- 
enced by the magnetic field in more detail, we depict the 
spatial dependence of (a) normalized local DOS for var- 
ious and magnitude of even and frequency pairings 
for (b) H/Ho=Om, (c) H/Ho=0.3, and (d) H/Ho=0.5 
dit E = y = in Fig. [TTl We see that at the core center, 
only odd frequency pairing exists while at the midpoint of 
the vortex lines, only even frequency pairing is present, 
which is reflected in the DOS as seen in Fig. [TT] (a). 
The magnitude of even and odd frequency pairings are 
reduced by increasing magnetic field. For higher energy, 
we also find a similar tendency (not shown). 



IV. SUMMARY 

In summary, we have studied pairing symmetry in the 
FFLO vortex and vortex lattice. We showed analytically 
that at the intersection point of FFLO nodal plane and 
vortex line, only even frequency pairing is present if the 
Zeeman splitting is negligibly small. With increasing Zee- 
man splitting, odd frequency pairing also emerges there. 
Therefore, the gap structure of the DOS at the intersec- 
tion point predicted in Refs. 13211331 can be regarded as a 
manifestation of the even frequency pairing. 

In the vortex lattice, at the core centers, only odd fre- 
quency pairing is present while at the midpoint of the 
vortex lines, only even frequency pairing appears irre- 
spective of the energy. Thus, the odd and even frequency 
pairings also form the lattice in the vortex lattice state. 

We have explained the electronic structure in vortex 
systems in terms of the pairing symmetry, even and odd 
frequency pairings, which has been understood with the 
quasiparticle picture to date. Our approach can be ex- 
tended to vortex system in other fields. Recently, vortex 
state has been realized in cold atoms. To study pairing 
symmetry in vortex state of cold atoms, as is done in this 
paper, would give an insight into vortex physics in cold 
atoms. 
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